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ABSTRACT 

The discovery of a low-mass star with such low metallicity as < 4.5 x 10~ 5 Zq reveals the critical 
role of dust in the formation of extremely metal-poor stars. In this paper we explore the effect of the 
growth of dust grains through accretion of gaseous refractory elements in very low-metallicity pre- 
stellar cores on the cloud fragmentation induced by the dust emission cooling. Employing a simple 
model of grain growth in a gravitationally collapsing gas, we show that Fe and Si grains can grow 
efficiently at hydrogen densities of ~ 10 10 -10 14 cm~ 3 in the clouds with metal abundances of —5 < 
[Fe, Si/H] < —3. The critical metal number abundances, above which the grain growth could induce 
the fragmentation of the gas clouds, are estimated to be ^4 cr it — 10 -9 -10 -8 , unless the initial grain 
radius is too large (> 1 /Ltm) or the sticking probability is too small (< 0.01). We find that even 
if the initial dust-to-gas mass ratio is well below the minimum value required for the dust-induced 
fragmentation, the grain growth increases the dust mass high enough to cause the gas fragmentation 
into sub-solar mass clumps. We suggest that as long as the critical metal abundance is satisfied, the 
grain growth could play an important role in the formation of low-mass stars with metallicity as low 
as 10~ 5 Z Q . 

Subject headings: dust, extinction - ISM: clouds - stars: formation - stars: low-mass - stars: Popu- 
lation II - supernovae: general 



1. INTRODUCTION 

Dust grains in the early universe are considered to be 
important agents to trigger the formation of low-mass 
stars in metal-poor environments (Omukai 2000; Schnei- 
der et al. 2003); the cooling of the gas through ther- 
mal emission of dust makes collapsing dense cores grav- 
itationally unstable, leading to the fragmentation into 
multiple sub-solar mass clumps at gas densities of 10 12 - 
10 14 cm" 3 (Tsuribe & Omukai 2006; Dopcke et al. 2011). 
This scenario has been recently supported by the discov- 
ery of a Galactic low-mass star, SDSS J102915+172927 
(Caffau et al. 2011), whose metal content is too low 
(Z < 4.5 x 10~ 5 Zq) to induce the fragmentation of 
star-forming clouds by metal-line cooling (see Klessen et 
al. 2012; Schneider et al. 2012b for details). 

The condition that realizes the dust-induced fragmen- 
tation depends on the amount of dust grains as well as 
their size distribution in pre-stellar clouds (Omukai et 
al. 2005; Schneider et al. 2006, 2012a). Schneider et al. 
(2012a) found that the formation condition of the low- 
mass fragments obtained by the numerical simulations 
is fully described in terms of the product of dust-to-gas 
mass ratio T> and geometrical cross section per unit dust 
mass S as follows; 

SV > 1.4 x 10" 3 cm 2 g" 1 ( -^-) 2 ( 19 ° H X h) 
6 V 10 K J V10 12 cm- 3 / VV 

where T gas is the temperature of the gas, and ch is the 
hydrogen number density. Treating self-consistently the 
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dust formation in the ejecta of supernovae (SNe) and the 
subsequent destruction of the dust by the reverse shocks, 
Schneider et al. (2012a) argued that the condition could 
not be satisfied in the collapsing star-forming clouds en- 
riched with metals and dust from the first SNe if a major- 
ity of grain formed in the SN ejecta are destroyed by the 
reverse shock (see also Schneider et al. 2012b). However, 
it could be possible that the accretion of gaseous refrac- 
tory elements released from dust grains in the shocked 
gas onto the surfaces of the SN dust surviving in star- 
forming clouds changes the mass and size distribution of 
the dust, and thus affects the thermal evolution of the 
collapsing cores. 3 

In this paper, we investigate the feasibility of grain 
growth in low-metallicity star-forming clouds to explore 
whether the grain growth can facilitate the formation 
of metal-poor low-mass stars. In Section 2, we describe 
the model of grain growth in collapsing dense clouds, 
and present the results of the calculations in Section 3. 
In Section 4, we estimate the critical metal abundances 
above which the grain growth could encourage the gas 
fragmentation into sub-solar mass clumps, and discuss 
the corresponding dust-to-gas mass ratio and total metal- 
licity. The conclusion is given in Section 5. 

2. MODEL OF GRAIN GROWTH IN METAL-POOR 
STAR-FORMING CLOUDS 

We consider the growth of dust grains in collapsing 
clouds that have been enriched with metals and dust 
grains produced by very early generation of SNe. Dust 
formation calculations by Nozawa et al. (2003) showed 

3 Hirashita & Omukai (2009) examined the coagulation of dust 
in collapsing clouds with a variety of metallicity. They found that 
the dust coagulation can proceed even at metallicity as low as 10 — 6 
Zq but does not have any impact on the thermal evolution of the 
star-forming clouds. 
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that various grain species condense in the unmixed ejccta 
of Population III SNe, and that Fe, Si, and C grains have 
relatively large average radii (> 0.01 /zm). Based on their 
dust models, Nozawa et al. (2007) investigated the evo- 
lution of dust in the hot gas swept up by the SN shocks 
and found that most of such large Fe, Si and C grains can 
survive the destruction by the reverse shock to be pre- 
dominantly injected into the early interstellar medium. 

Being motivated by these studies, and to simplify the 
chemistry, we consider that Fe and Si grains composed 
of only one element grow through accretion of Fe and 
Si atoms in the gas phase, respectively. 4 We assume 
that dust grains are spheres, all of which have a single 
initial radius, although they might be expected to have 
the distribution of grain sizes. Suppose that the number 
density of a given refractory element i at a time to is 
c ;,o = Cj(io) = AiCnfl, where At is the number abundance 
of the element i relative to hydrogen. In order to specify 
the number abundance of pre-existing seed grains, we 
introduce a parameter /^o (>0), which is defined as the 
number fraction of the elements i originally locked in 
dust grains (i.e., condensation efficiency at to). Then, 
the number density of the dust whose initial radius is 
r ij0 is described as nf^ st = fi,oCi,o(ai,o/n,o) 3 , where a ij0 
is the hypothetical radius of an atom in the dust phase. 

The time evolution of the number density Ci(t) of an 
element i in the gas clouds collapsing with the timescale 
of free-fall is given by 

Ci(t) = c ifi ^1 - -^-^j , (2) 

where t is the elapsed time from to, Tq = 
(37t/32G/itohCh,o) 1 ^ 2 is the free-fall time at the density 
ch.o with the gravitational constant G, the mean molec- 
ular weight fi, and the mass of a hydrogen atom ran- 
Once the grain growth activates, the gaseous atoms are 
consumed, and the number density cf as (£) of an element 
i in the gas phase at the time t can be written as 

cr(t) = c l (t){i-f h0 [n(t)/r lfi f}, (3) 

with ri(t) being the radius of the i-th grain species at t. 
Equation (3) is reduced to 

f i (t) = l-Y i (t)=f ifi Xf(t), (4) 

where fi(t) is the condensation efficiency at t, Yj(t) = 
cf as (t)/a(t) represents the depletion of the gaseous 
atoms due to grain growth, and Xi(t) = ti (t) /?'.;, o- 

In the dense clouds where almost all gaseous atoms are 
neutral, the growth rate of grain radius is given by 

where Sj is the sticking probability of the gaseous element 
i incident onto grains, k is the Boltzmann constant, and 
rrii is the mass of the clement i. The gas temperature T gas 
is assumed to be constant during the evolution of clouds 

4 The growth of C grains may not be expected in dense clouds. 
This is because at high gas densities considered in this paper (ch > 
10 8 cm -3 ), all C atoms can be bounded in CO molecules in gas 
clouds with metallicity scaled by the solar abundance (Omukai et 
al. 2005). 



TABLE 1 

Numerical Values Used in the Calculations 



Numerical Values 


Explanation of Symbols 


s — 1 


■ i ■ u un- 
sticking probability 


T gas = 10 3 K 


gas temperature 


H = 2.18 


mean molecular weight 


a FCi0 = 1-441 A 


radius of a Fe atom in the solid phasc a 


a Si! o = 1.684 A 


radius of a Si atom in the solid phase a 


mp c = 56jtih 


mass of a Fe atom 


msi = 28rrtH 


mass of a Si atom 


A Fe , = 3.26 x 10~ 5 


solar abundance of Fe relative to H b 


A Si ' e = 3.58 x 10~ 5 


solar abundance of Si relative to H b 


PFe = 7.90 g cm" 3 


bulk density of Fe c 


PSI = 2.32 g cm" 3 


bulk density of Si c 



a Nozawa et al. (2003) 

b Anders & Grevesse (1989) 

c Nozawa et al. (2006) 



in this study. The supcrsaturation ratio Si is a function 
of the dust temperature Tdust- Since Td us t/T gas <C 1 and 
Si ^S> 1 under the condition considered here (Dopcke et 
al. 2011), Equation (5) is reduced to 

by introducing ( T i S o) _1 = 

Si47ra? (fcr gas /27rm i ) 1 ' 2 AiCH,o/3r i;0 . Then, inte- 
gration of Equation (6) leads to the ratio of grain radius 
to the initial one Xi(t); 

2r n ff f u Y(u') 
Xi(t) = l + -± / JT^du', (7) 

T ifi J (1 - u') 

with u = 4/2tq . 

In principle, by solving Equations (4) and (7) for a 
given set of /^o, ri t o, and Aj, we can calculate the time 
evolution of fi(t) and r^(t) = ri^Xi(t). The values of 
the other parameters necessary for the calculations are 
summarized in Table 1. As is shown later, for the metal 
abundances considered in this paper, the grain growth 
operates at high gas densities of ch > 10 10 cm -3 , where 
the gas temperature is expected to be in the range of 
500-2000 K (Dopcke et al. 2011). Thus, the calculations 
are started from ch,o = 10 8 cm" 3 , with T gas = 10 3 K. 
Note that the results of calculations are not sensitive to 
Tgas as long as the above range of T gas is considered. 

3. RESULTS OF CALCULATIONS OF GRAIN GROWTH 

Figure 1 depicts the growth of Fe grains with /f c ,o = 
0.1 and rF e ,o = 0.01 /jm; the time evolutions of the con- 
densation efficiency f-Feif) and grain radius r-p e (t) versus 
hydrogen number density cn(t) for [Fc/H] = —5, —A, and 
-3, which correspond to A Fc = 3.26 x 10~ 10 , 3.26 x 10~ 9 , 
and 3.26 x 10 -8 , respectively, with the solar abundance 
by Anders & Grevesse (1989). We can sec that the 
grain growth activates efficiently even in the gas clouds 
with [Fe/H] = —5, and the grain radius finally reaches 

a constant value rF e ,o(l//Fe,o) 1 ' 3 by consuming up all 
gaseous Fe atoms. However, the gas density at which a 
considerable fraction (/p c ~ 0.5) of Fe atoms is locked 
up in dust grains is higher for a lower Fe abundance; 
c H ^ 10 10 , 10 12 , and 10 14 cnL 3 for [Fc/H] = -3, -4, 
and —5, respectively. Also, the gas density at which /p e 
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Fig. 1. — Time evolutions of the condensation efficiency /F e (t) 
(red) and grain radius rp c (t) (blue) through the growth of Fe grains 
with the initial dust abundance /p e q = 0.1 and the initial grain 
radius rp e rj = 0.01 (im as a function of hydrogen number density 
cjj(t). Dot-dashed, solid, and dashed lines depict the results for 
the collapsing clouds with the Fe abundances of [Fe/H] = —3, —4, 
and —5, respectively. 

reaches ^0.5 is two orders of magnitude higher (lower) 
for rpo.o = 0.1 (0.001) /im than that for rp c .o = 0.01 
/im, though not presented in the figure. These behaviors 
of grain growth in the collapsing gas clouds can be seen 
from Equation (9) (see below); for given values of /^o 
and /;,*, cn(t) oc (r ij0 Mi) 2 if Ca{t)/ca,Q > 1- Hence, 
for a fixed fi,o, one order of magnitude higher r^o or one 
order of magnitude lower Ai is compensated with two 
orders of magnitude higher cjj- 

Figure 2 shows the time evolutions of /si(t) and r$i(t) 
of Si grains with rs;.o = 0.01 /im for different initial 
dust abundances of /si,o = 0.1, 0.01, and 0.001. Here, 
the total abundance of Si atoms is set to be [Si/H] = 
-4 (A Si = 3.58 x 10" 9 ). For f ifi = 0.1, the growth 
of Si grains proceeds somewhat earlier than that of Fe 
grains for [Fe/H] = —4, and the condensation efficiency 
increases to fsi = 0.5 at ch — 10 11 cm -3 . In the cases of 
lower /si,o, higher gas densities are needed for achieving 
some level of the condensation efficiency, although the 
final grain radii are larger when all Si atoms are tied up 
in dust grains. 

4. DISCUSSION: CRITICAL METAL ABUNDANCES 

As shown in the last section, the grain growth can 
activate even in metal-poor star-forming clouds whose 
metallicity is only [Fe, Si/H] ~ —5. However, in or- 
der that the grain growth affects the thermal evolution 
of collapsing cores, it must become effective before the 
cloud density increases to ch = 10 12 -10 14 cm -3 , where 
the optical depth becomes high enough to suppress dust 
emission cooling, and as a result the gas fragmentation 
is expected to occur (e.g., Schneider et al. 2012a). 

The metal abundance above which the grain growth 
becomes important can be estimated by requiring that a 
certain fraction fi » of the element i should be locked up 
in dust grains at a given hydrogen number density ch.*- 
Noting that F;(t) = 1 - f ifi Xf(t) = 1 - fi(t), Equation 
(7) can be rewritten as 



Ki(fi,0, fi, 



X, 



dX, 



c H (t) (cm ) 

Fig. 2. — Time evolutions of the condensation efficiency /si(i) 
(red) and grain radius rgi(t) (blue) through the growth of Si grains 
°f r Si = 0-01 A™ 1 m the collapsing clouds with the Si abundance 
of [Si/H] = —4. Solid, dashed, and dot-dashed lines represent the 
results for the initial dust abundances of /si o = 0.1, 0.01, and 
0.001, respectively. 
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du' 



2r ff 
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(8) 



where = (/ 4 ,*//i, ) 1/3 and u* = 1 - (c h ,o/ch,*) 1/2 - 
From this equation, we can derive the critical metal 
abundances above which the grain growth could facili- 
tate the fragmentation of the clouds as follows 



Aj 



(1.0 - 2.5) x 10~ 9 Ki 



0.01 /im 



10 12 cm" 3 
ch,* 



1 - fifiXi 



using ch,o/ch,* -C 1. The numerical factor 2.5 (1.0) 
corresponds to Fe (Si) grains. The function Ki in- 
creases with increasing /j * and/or decreasing /j q; for 
0.2 < ft* < 0.8, Ki = 0.3-2.4 (5.2-15) at f i0 = 0.1 
(fi,o = 0.001). 

Figure 3 presents the critical abundances of Fe and Si 
in the form of [X/H] versus fi$ for r^o = 0.01 /im; the 
dot-dashed, solid, and dashed lines give the abundances 
necessary for the condensation efficiency to reach 0.8, 
0.5, and 0.25, respectively, at ch,* = 10 12 cm' 3 . As ex- 
pected, higher metal abundances are needed for attain- 
ing higher and/or for lower f^g. For the case of 
fi,* = 0.5, the critical abundances of Fe and Si spans the 
ranges of -4.11 < [Fe/H] < -3.19 and -4.54 < [Si/H] 
< —3.62, respectively, for the range of 0.1 > /).o > 0.001. 
It would be interesting to mention that the above range 
of [Si/H] covers the abundance of Si observed for SDSS 
J102915+172927 ([Si/H] = -4.27, Caffau et al. 2011). 
This could suggest that the growth of Si grains might 
have worked in the parent cloud of this star. 

Here we present how the fragmentation condition in- 
duced by grain growth depends on the unknown param- 
eters such as ri : o and Sj. The time duration At s for 
which the gas density increases from ch,o to ch,* by free- 
fall is given by At s = 2t§ (1 - x c H .. m.- :■• On the 
other hand, the time duration At gs for which the con- 
densation efficiency /^o at ch,o increases up to at 
ch,* through grain growth is derived from Equation (8) 
as (Am)- 1 = (A' i T gs )- 1 + (2T ff )- 1 . Since At s /Am > 1 
is required for the fragmentation, we can obtain the con- 
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Fig. 3. — Critical abundances of Fe (upper panel) and Si (lower 
panel), for which the grain growth can lock up a fraction = 
0.8 (dot-dashed), 0.5 (solid), and 0.25 (dashed) of Fe and Si atoms 
in dust grains at the hydrogen number density ch,* = 10 12 cm -3 . 
The horizontal axis is the initial dust abundance /^Oi an d the initial 
grain radius is set to be r^o = 0.01 (mi. 

dition 



(1-0 - 2-5) 



A, 



2.5 x 10- 9 
ch,* 



0.01 urn 



1.0 



>1.0, (10) 



where the numerical factor 1.0 (2.5) is for Fe (Si) grains. 
This inequality demonstrates that the fragmentation 
condition is achieved more easily for the initial grain radii 
smaller than 0.01 /im. In contrast, if r^o > 1.0 /im or 
Si < 0.01, the grain growth can no longer induce the 
fragmentation for A, < 2.5 x 10" 7 ([Fe, Si/H] < -2). 
We also note that a smaller /^o producing a larger Ki 
acts against the gas fragmentation by grain growth. 

Next we consider the dust-to-gas mass ratio to sec 
whether the condition for the dust-induced fragmenta- 
tion given in Equation (1) can be met. In the con- 
text of this paper, we suppose that the product of S 
and T> without grain growth is given as (<SZ?)i iC rit = 
3/i,o^4i.crit^i/4piri ! oA tm H, adopting the critical metal 
abundance evaluated in Equation (9). On the other 
hand, the product resulting from the grain growth is 
given by {SV)^ = (SV)^ clit (f h */ f h0 ) 2 / 3 . Figure 4 
shows the dependence of (<SX>)i, C rit and (<S2?)i,* on /^o, 
adopting = 0.5. We can see that (<SZ?)i jC rit is we H 
below the minimum value required for the dust-induced 
fragmentation (dot-dashed line in Fig. 4), whereas 
(SV)i^ exceeds this value. This indicates that, even if 
the destruction by the SN reverse shock results in a lower 
ST> than the criterion for the dust-induced fragmenta- 
tion, the grain growth can enhance ST> in the clouds 
and can enable the gas fragmentation into sub-solar mass 
clumps. Note that the results in Figure 4 are indepen- 
dent of the initial grain radius r^o since Si oc and 
T>i oc A iiCI i t oc rift. 

Finally, we relate the critical abundance to the to- 
tal metallicity Z. Equation (9) suggests that the crit- 
ical metal abundances are generally in the range of 
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Fig. 4. — Products of dust-to-gas mass ratio T> and geometrcal 
cross section per unit dust mass S for Fe (blue) and Si (red) as 
a function of the initial dust abundance /^j. The dashed lines 
depict (ST>)i cr ; t corresponding to the critical abundances given 
by the solid lines in Figure 3, whereas the solid lines depict the 
resulting dust-to-gas mass ratio after the grain growth (SZ')i ». 
The horizontal dot-dashed line indicates the minimum value above 
which the dust emission cooling causes the gas fragmentation into 
low-mass clumps (Schneider et al. 2012a). 

A crit ^ 10" 9 -10~ 8 , depending on / ij0 , /,,*, r i)0 , and grain 
species. By representing the mass ratio of refractory el- 
ements condensible into dust grains to the total heavy 
elements as 7?., the metallicity Z can be related to j4 cr i t 
as 1ZZ ~ A CI i t (Hrei I n) with /i rc f being the mean atomic 
mass of refractory elements. Then, we have 



(5 - 50) x 10" 



02 



(11) 



where we use /x re f//i = 20 and Zq = 0.02. Equation (11) 
implies that, if the grain growth does work efficiently, it 
can drive the gas fragmentation of low-mass clumps in 
the star-forming clouds enriched with metallicity ~10 -5 
Zq. In other words, as long as the abundance of a 
given refractory element satisfies the critical abundance 
in Equation (9), the formation of hyper- metal-poor low- 
mass stars with the metallicity lower than Z ~ 4.5 x 10~ 5 
Z & observed in SDSS J102915+172927 could be possible. 

5. CONCLUDING REMARKS 

We have investigated the growth of dust grains in 
metal-poor proto-stellar clouds. Our simple model shows 
that the grain growth can operate efficiently even in col- 
lapsing dense cores with metal abundances as low as 
[Fe, Si/H] ~ —5. We also present the critical metal 
abundances above which the grain growth could affect 
the fragmentation process of collapsing gas clouds. This 
abundance is estimated to be A cr - lt — 10 -9 -10 -8 , which 
suggests that the formation of low-mass stars with metal- 
licity of ~10 -5 Zq can be possible. We conclude that 
even if the initial dust-to-gas mass ratio does not satisfy 
the condition required for the dust-induced fragmenta- 
tion, the grain growth can increase the dust-to-gas mass 
ratio high enough to facilitate the formation of metal- 
poor low-mass stars. 

Our results suggest that if grain growth is considered, 
the formation of low-mass protostars can occur not only 
at very low metallicity but also at higher metallicity. 
The final mass of a newly born star is determined by 
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the accretion of the surrounding gas onto the protostars 
(McKee & Ostriker 2007 and references therein), and its 
mass accretion rate would be regulated by the mass of 
the central protostar induced by the grain growth. Thus, 
the grain growth in collapsing clouds might be a funda- 
mental physical process to control the stellar initial mass 
function in the present universe. 

It should be mentioned that we have considered only 
the growth of single-component Fe and Si grains with a 
single initial radius. However, it might be possible that Si 
and Fe atoms condense as silicates or oxides in an oxygen- 
rich gas. Since the growth of such compound grains with 
no monomer molecule has been usually treated by consid- 
ering Si or Fe element as a key element (e.g., Zhukovska et 
al. 2008), the mass and radius of dust given in this paper 
are considered to be lower limits. On the other hand, 
the timescale of grain growth is sensitive to the initial 
grain radius (Hirashita & Kuo 2011). Thus, the effect of 
the initial size distribution as well as the growth of corn- 
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pound grains should be explored. Furthermore, we have 
assumed the sticking probability of Sj = 1 and a con- 
stant gas temperature T gas = 10 3 K during the collapse 
of the clouds. In particular, too low sticking probabilities 
(si < 0.01) may prevent the grain growth from becoming 
efficient for the metal abundances of A% < 2.5 x 10~ 7 . We 
note that our conclusions obtained by the simple model 
should be confirmed by more sophisticated simulations 
of the thermal evolution of star-forming clouds involving 
grain growth. 
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